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Abstract 

The PYTHIA program can be used to generate high-energy-physics events with sets of 
outgoing particles produced in the interactions between two incoming particles. The ob- 
jective is to provide a representation, as accurate as possible, of event properties in a 
wide range of reactions. One such reaction, that is not well understood is Diffraction. 
Among the several alternative approaches that have been proposed, in PYTHIA, we fol- 
low a fairly conventional Pomeron based one, but fully integrated to use the standard 
PYTHIA machinery for multiple interactions, parton showers and hadronization. This 
note reports the development in PYTHIA in the way diffraction is modeled without pro- 
viding specific details for usage. Results are compared with an alternative event generator 
called PHOJET. The code and further information may be found on the Pythia web page: 



http : //home . thep . lu . se/~torbj orn/Pythia. html 
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1 Diffraction and the Pomeron 

Hadronic processes can be classified as being either soft or liard. Soft processes, that domi- 
nate hadronic scattering cross-sections, are characterised by an energy scale of the order of the 
hadron size {lim ^ 200 MeV) pp. The hard sector is described very well by Perturbative QCD 
(pQCD) scatterings. However, pQCD is inadequate to describe soft processes, as a large scale 
makes the coupling constant (a^) large enough to make the higher order terms non- negligible, 
thus making the process intrinsically non-perturbative. 

I. Pomeranchuk predicted that if the total cross-section behaves asymptotically like a power 
of Ins, then the particle and anti-particle cross sections become asymptotically equal [2]. The 
exchange of a Regge trajectory that ensures this behaviour was first introduced by Gribov [5]. 
The particles on this trajectory are virtual and have the same internal quantum numbers as 
the vacuum. The effective summation of particles on this trajectory is known as the Pomeron 
(P). In QCD, the Pomeron is regarded as a colourless and flavourless multiple gluon [5] or a 
glueball exchange. 

1.1 Classification 

In proton-proton (pp) (or more generally hadron-hadron) scattering, interactions are classified 
by the characteristics of the final states. Interactions can either be elastic or inelastic. In elastic 
scattering {pi +P2 ^ Pi +P2)) both protons emerge intact and no other particles are produced 
as shown by the pink dots in figure [TJ The LHC cross-section (at ^/s = 14 TeV) for elastic 
scattering is ~30mb [5]. 

Colliding hadrons are colour singlets. As they approach each other, they may exchange a colour 
octet gluon, making each hadronic cluster a colour octet. As they move apart, they need colour 
lines connecting them. To be able to separate into two separate systems, they need to exchange 
another gluon and become colourless. However, the final state need not be identical to the ini- 
tial state. Such processes are called inelastic. When colour octets move apart, colour lines are 
stretched between them. Given time, this system gets complex and multi-particle production 
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Figure 1: Diagram for elastic scattering and cj) vs t] pfot showing the distribution of products after the 
interaction. 



occurs. 



Elastic scattering can be achieved via the exchange of a glueball-hke Pomeron. In elastic scat- 
tering, the final state and initial state particles are identical. The exchange of gluons can excite 
a hadron. This can result in the outgoing state preserving the internal quantum numbers of 
the incoming particles but having a higher mass. This is known as quasi-elastic scattering. 

Inelastic collisions can be diffractive. There are several possible descriptions of diffraction, al- 
lowing several alternative approaches. The approach discussed here is one described by Regge 
theory [B] in terms of the exchange of a Pomeron. One of the alternative approaches which does 
not use the concept of a Pomeron or Regge phenomenology is called the soft colour interaction 
model. It is described by Ingelman in [7]. 



A diffractive reaction is one in which no internal quantum numbers are exchanged between the 
colliding particles. Diffraction occurs when the exchanged Pomeron interacts with the proton 
to produce a system of particles referred to as the diffractive system. In diffractive scattering, 
the energy transfer between the two interacting protons remains small, but one or both pro- 
tons dissociate into multi-particle final states with the same internal quantum numbers of the 
colliding protons. 

If only one of the protons dissociates then the interaction is Single Diffractive (SD) {pi +P2 -^ 
p[ + X2 or pi + P2 -^ Xi + P2). The dissociated proton is shown as a spray of blue dots (par- 



tides) and the non-dissociated proton as the pink dot in figure [2l The LHC cross-section (at 
yi = 14TeV) for SD is ~ 10 mb [5j. 




Figure 2: SD diagram and a window showing a rapidity gap between —10 < rj < 3.5. 

If both the colhding protons dissociate, then it is Double Diffractive (DD) {pi +p2 —> Xi +X2) 
as seen in figure El The LHC cross-section (at ^/s = 14TeV) for DD is ~ Tmb |5]. 
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Figure 3: DD diagram and window showing a rapidity gap between —3.5 < 77 < 4. 

A different topology becomes possible with two Pomerons exchanged, namely Central Diffrac- 
tion (CD) {pi + P2 -^ p'l + X + P2) or Double Pomeron Exchange. In this process, both the 
protons are intact and are seen in the final state (as two pink dots seen in figure S]). The LHC 
cross-section for CD is ~ 1 mb [5] . 



In Non-Diffractive (ND) interactions there is an exchange of colour charge and subsequently 
more hadrons are produced. This is shown in figure |5l ND interactions are the dominant 
process in pp interactions and are expected to be ~60% of all interactions at the LHC with a 
cross-section of ~65mb (at y/s = 14TeV) [5]. 




Figure 4: CD diagram and window showing two rapidity gaps between — 10 < 77 < —2.5 and 2.5 < 77 < 10. 





Figure 5: The diagram for an ND process. The rapidity window on the right shows that there is no rapidity 
gap. 

A consequence of the Pomeron hypothesis is that the cross-sections of pp and pp diffractive 
scattering should be equal at high enough energies [S]. This is because the Pomeron has the 
quantum numbers of the vacuum, so its couplings to the proton and anti-proton are equal. 
The total pp cross-section is given by equation [1] where "misc" here is CD and multiple Pomeron 
exchange. The cross-section for multiple Pomeron exchange is ^ 1 mb [5] . 



0"tot = CTel + CTinel = (^e\ + CTdiS + CT-^D = (^e\ + OSD + CTDD + (^n 
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1.2 Kinematics 

In a QCD approach, a partonic description of a Pomeron, as described in [9j is commenly used. 
Distributions of partons in particles are characterised by Parton Distribution Functions (PDF). 
A PDF fi{x, Q"^) gives the probability of finding a parton i with a fraction x of the momentum 
of the parent beam particle, when probed at a scale of Q^ . PDFs are parameterisations of 
experimental data. Diffractive hard scattering is used to resolve the partonic structure of the 
Pomeron [TU]. 



Different alternative factorizations of tlie partonic structure of the Ponieron exist. A model for 
diffractive hard scattering is described in [TT]. In this type of factorisation, firstly a Ponieron is 
emitted from a proton pi (at the upper vertex in figure [6]) with a momentum transfer squared 
given by 



t = iPi- p'i) 
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Figure 6: Exchange diagram for single diffraction. 

Then this emitted Ponieron interacts with the other proton, pj at the lower vertex, with a 
transfer of momentum between constituent partons. The system X that is produced in this 
interaction is called the diffractive system. There is a large rapidity gap (LRG) between the 
out-going proton and diffractive system X. This introduces the concept of a Pomeron flux in 
a proton {p) /p/p(a;p,t) (in this case /p/pi(xp, t)), where xp is the fraction of the proton's 
momentum carried by the emitted Pomeron, and diffractive PDFs (DPDF). The Pomeron flux 
describes the probability that a Pomeron with a given value of xp and t couples to the proton. 
In the massless limit, 

xp = E^/Ep (3) 



where E-p and Ep are the energy of the Pomeron and the proton to which it was coupled to 



respectively. The fraction of the Pomeron's momentum carried by its constituent gluon (g) (or 
quark (q)) is given by 

Xg (or q) = Eg (or q)/-E'p (4) 

where E is the energy of the gluon (or quark) . 

The diffractive hard pp scattering cross-section can be written as 

dajpp^p + X) ^ ^^^^^^^^ ^^ daipF^X) ^ ^^^ 

dxpdtdxidx2dt ^ ^ ■' dxidx2dt 

Pfiux 

The second term in equation O is the proton-Pomeron hard scattering differential cross-section. 
It can be assumed to factorize as in equation [G) 

daipV ^ X) ^2\ r / ^2\d'd' ,„. 

, , ,2 = fpi/pi^i^ Q )fp2/Ax2, Q )-j- (6) 

dxidx2dt dt 

Here, fpi/p{xi,Q^) and fp2/p{x2,Q'^) are the proton and Pomeron PDFs with partons pi and 
P2 having momentum fractions Xi and X2 of the proton and Pomeron respectively, da /dt is 
the corresponding hard scattering cross-section for that subprocess. Because of the inherent 
non-perturbative effect in a QCD binding state, PDFs (and DPDFs) cannot be obtained by 
perturbative QCD from first principles. The known PDFs (and DPDFs) are instead obtained 
by using fits to experimental data. The DPDFs used here are obtained using the DGLAP 
evolution equations [12]. The invariant mass of the diffractive system X, also known as the 
diffractive mass, is given in terms of the overall collision cms energy y/s by 

Ml = xps. (7) 

Experimentally, diffractive reactions are characterised by a large (non exponentially-suppressed) 



rapidity gap in the forward region, e.g, xp < 0.1. In other words, there is a large separation in 
rapidity between the quasi-elastically scattered proton and the diffractive system, in which no 
particles are detected. A few ND events may also display a large rapidity gap due to multiplicity 
fluctuations but their number is exponentially suppressed with increasing rapidity gap. 
Another exchange mediator called the Reggeon R [13] is needed to reproduce experimental 
data of diffractive processes and total cross-sections successfully. Thus, there are two exchange 
mediators: Reggeons and Pomerons. Reggeon exchange fits data at relatively lower energies 
(high xr), as Reggeons couple to valence quarks of a proton, which carry a large fraction of the 
Pomeron's momentum x. At high energies, the incoming protons "pass by" so quickly that it is 
mainly the sea quarks that interact. On the other hand, Pomeron exchange fits the data only 
at higher energies (low xp), as a Pomeron couples to gluons (and sea quarks). Already at ISR 
energies {^/s = 63[GeV]), the exchange mediator was predominently the Pomeron. Thus, the 
higher the collision energy, the more important is the role of the Pomeron. The sum of these 
two trajectories describe the total pp cross-section. 

Based on these theories about the Pomeron, a model for diffraction has been constructed and 
implemented in Pythia [l^, resulting in a complete final state. 



2 Diffraction in PYTHIA 
2.1 PYTHIA 6 

The development of series 6 of PYTHIA written in Fortran 77 began in 1997. Although 
there was significant development from one version to the next, the description of diffraction 
remained the same (in the two versions 6.2 and 6.4). In this section a description of the 
diffractive processes in PYTHIA 6 is presented. Also presented is a comparison of diffractive 
kinematic distributions produced using PYTHIA and an alternative Monte Carlo generator 
called PHOJET 



2.1.1 Event Generation 

The total hadronic cross-section for AB — )■ anything, a^ is calculated using the Donnachie 
and Landshoff parameterization [8]. In this approach, the total cross-section appears as a sum 
of a Pomeron term and a Reggeon term, as seen in equation [HI 



t^{s)=X^^s^ + Y^^s-^ (8) 



cr, 



The powers e for the Pomeron term and rj for the Reggeon term are expected to be universal, 
while the coefficients X and Y are specific to each initial state. Those parameterizations not 
provided in [8] have been calculated in [16j. Cross-sections for elastic, single and double diffrac- 
tive events are provided, but higher diffractive topologies like central diffraction are neglected. 
The diffractive cross-sections and event characteristics are described by a model by Schuler and 
Sjostrand found in [161 [TTj. The non diffractive cross section is given by "whatever is left". 

In the Schuler-Sjostrand model, when the square of the momentum transfer t is not too large, the 
differential elastic cross section can be approximated by a simple exponential fall-off with respect 
to t. Diffractive cross-sections have an inverse diffractive mass squared (M^) dependence and 
an exponential dependence on t. The simple dM'^/M'^ form is modified by the mass dependence 
of the slope of the t distribution (with co-efficients -Bsd(XB), -Bsd(AX) and Bdd) [ISJ. These Regge 
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formulae for diffraction are supposed to hold in certain asymptotic regions of the full phase 
space. For example, in p + p — ;• p + M, 1^^^^^ ^ mp(M^ — ni^/ s [17J. Having a lower cut 
on t of the order of ra^ implies M^ — m^ < 0.15s. But there will be diffraction even outside 
these regions where the Regge formulae were derived. Due to the lack of a theory that predicts 
differential cross-sections at arbitrary t and M^ values, the above Regge formulae are used 
everywhere along with "fudge factors" F^.^ and Fjd in equations [3 and (TDl The form of these 
factors is given by equations 7.76 of [18] to give sensible behaviour in full phase space. These 
factors suppress production close to the kinematical limit and in the case of double diffraction, 
also suppresses configurations where the two diffractive systems overlap in rapidity space. These 
"fudge factors" also give a broad enhancement in the production rate in the resonance region 
up to about 2 GeV. This gives a smeared-out version of exclusive states, rather than listing 
them all out individually. 
Diffractive cross-sections are given by equations O and [TUl 

didM- = 16^^^^^^^ M^ exp(5,d(XB)t)F3d (9) 

dcTddi-s) _ 93F a n 11 



dtdMldMl 167r '^^^'^^^ Ml M'. 



/3AP/3BP7T2 7T2eXp(5ddt)Fdd (10) 



The couplings /3 are related to the Pomeron term of equation [HI The triple Pomeron coupling 
(73P is determined from single diffractive data. The diffractive mass spectrum M ranges from 
0.28 GeV (~ 2r?7,^) above the mass of the diffracted hadron, to the kinematic limit. The ex- 
ponential slope parameters -Bsd or dd are assumed to have a logarithmic dependence on 1/M^. 
The kinematic range of t depends on the masses of all incoming and outgoing systems involved. 
More information and the equations can be found in section 7.7.1 of the PYTHIA 6.4 manual 

m. 



Diffractive cross-sections have been integrated for a set of centre of mass (CM) energies, start- 
ing at 10 GeV. The results have been parameterized in section 4 of [ITj. Once the process is 
selected using this parameterization, M and t are generated using equations [9] and [TOl 



2.1.2 Particle Production 

Once the process is selected and the kinematic variables are determined, the products of the 
collision are generated. The handling of this production depends on the relative value of the 
diffractive mass M. If M < 1 GeV above the mass of the incoming particle, the system is 
allowed to decay isotropically to a two-body system. For a more massive diffractive state, the 
system is treated as a string with the quantum numbers of the original hadron. Two alternative 
ways of stretching the string are considered. 

There is both a gluonic and a quark contribution. When an incoming hadron is diffractively 
excited, either a valence quark or a gluon is "kicked out" of it. If the Pomeron couples to 
a valence quark of the non-diffracted proton, the string (the pink dashed lines in figure [7]) is 
stretched between the struck quark and the remnant diquark (or antiquark) of the diffractive 
system, seen in figure [Tl^a). This configuration dominates at small M. The alternative is when 
the interaction is with a gluon of the non-diffracted proton. The string is stretched from a quark 
in the diffractive state to a gluon, and then back to a diquark (or antiquark). This gives rise 
to a hair-pin structure as seen in figure [7[b). In PYTHIA 6 the ratio of the two contributions 
can be set. 



^ ^mQQQQQQSh 



(a) (b) 

Figure 7: String being stretched in diffractive processes - (a) P couples to a valence quark and (b) P couples 
to a gluon 
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2.1.3 PYTHIA 6.214 vs PHOJET 1.12 

A study comparing the pseudorapidity (77), charged particle density (dNch/drj) and transverse 
momentum (px) distributions in PYTHIA 6.214 and PHOJET 1.12 at CM energy 7TeV is 
shown below. ND and SD spectra are compared to analyse the difference in the diffractive 
part. 
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Figure 8: 77 distributions for (a) ND and (b) SD events at 7TeV comparing PYTHIA6 and PHOJET. 



A comparison of figures Eia) with [H](b) and M^a) with M^^) shows that although the multiplicity 
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Figure 9: Multiplicity distributions for (a) ND and (b) SD events at 7TeV comparing PYTHIA6 and PHOJET. 
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Figure 10: pr distributions for (a) ND and (b) SD events at 7TeV comparing PYTHIA6 and PHOJET. 
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spectra for ND events in PYTHIA and PHO JET are similar, high multiphcity SD events are not 
generated by PYTHIA. Similarly, the pt spectra in figures [TUlfa) and fTUT b) show that PYTHIA 
lacks a hard diffractive part. 

2.2 PYTHIA 8 before 8.130 

PYTHIA 6.4 was the last version of PYTHIA to be coded in Fortran 77, followed by a switch 
to C++ with version 8.1. The mechanism for diffractive scattering works in almost the same 
way as in PYTHIA 6. The only difference lies in the particle production. In PYTHIA 8.1 the 
ratio of the probability of the Pomeron coupling to a quark (-P(q)) and the Pomeron coupling 
to a gluon (-P(g)) is given by equation [11] A^ in this equation is a normalization factor and p 
(default value = 1) is a user-defined power. This introduces a mass dependence on the ratio of 
the two couplings, enabling the gluonic contribution to dominate at higher diffractive masses. 

^ = A (11) 

P(g) MP ^ ' 

2.3 PYTHIA 8.130 

In the versions of PYTHIA following PYTHIA 8.130, diffraction is modelled based on the 
Pomeron approach described in section 1.2 . Pomeron-proton collisions are modeled at a re- 
duced CM energy, then fully integrated in such a way that it uses the standard PYTHIA 
machinery for multiple interactions, parton showers and hadronization. This is the approach 
pioneered in the POMPYT program [TU] . 

2.3.1 Event Generation 

Diffractive cross sections are determined in exactly the same way as described in section 12.1.11 
However, in addition to the Schuler-Sjostrand model for picking M and t, three other parame- 
terizations of the Pomeron flux have been implemented. 
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1. Bruni and Ingelman [20]: based on a critical Ponieron giving a mass spectrum close to 
dM'^ jM'^ . The t distribution is the sum of two exponentials. 

2. Berger et al. [21] and Streng [22]: a conventional Pomeron description but with values (from 
the RAPGAP manual [2H]) updated to a super-critical Pomeron. This gives a stronger peaking 
towards low-mass diffractive states. The t slope is still exponential and depends on M. 

3. Donnachie and Landshoff [23]: similar to the Streng-Berger parameterization, but with a 
power law distribution for t. 

A comparison of the 4 different Pomeron fluxes are seen in flgure [Til 
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Figure 11: Differential cross-sections using different Pomeron fluxes. 

The Schuler-Sjostrand model is currently the only one which provides a separate t spectrum 
for DD. 
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2.3.2 Particle Production 

PYTHIA by default only allows collisions with CM energy above 10 GeV. But the diffractive 
mass spectra extend down to about 1.2 GeV. A perturbative description at this scale is not 
possible, giving rise to a separate handling of low and high masses. For M < 10 GeV, the non- 
perturbative description with longitudinally stretched strings, as described in sections |2 . 1 .21 and 



12.21 is implemented. 



High-Mass Diffraction In the mass range 10 GeV < M < ^/s, a perturbative description 
is implemented. The probability for this description is given by 

-Ppert = 1 - exp((mdiflf - mi„in)/"^width) 

where mmin and m„idth are free parameters. The default value of mmin is set at 10 GeV so that 
-Ppert vanishes when M < 10 GeV. 

The perturbative description involves using PDFs for the Pomeron that are not well known. 
PYTHIA 8.130 provides a selection of five PDF sets. 

1. Q^-independent parameterizations of the form given by equation [T2j Here A^ is a normal- 
ization factor that ensures unit momentum sum and a and h can be different for the quark and 
gluonic content of the Pomeron. In this PDF set, the momentum fraction of gluons and quarks 
can be freely mixed. Additionally, the production of s quarks can be suppressed relative to u 
and d quarks, with quarks and anti-quarks being equally likely to be produced. 

xf{x) = N^^,x''{l-xf (12) 

2. A Pomeron can be described by the PDF for a pion. A few PDF sets exist, one of which 
is built into PYTHIA. The others can be accessed from LHAPDF [2S]- Parameterizations are 
given for Tr"*"; 7r~ is obtained by charge conjugation and vr'' by averaging. 
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3. The HI 2006 Fit A parameterization is a Q^-dependent set. This is based on a tune (tune 
A) to HI data on inclusive diffractive cross-section, described in section 5.3 of 



4. The HI 2006 Fit B parameterization is another Q^-dependent set based on tune B to the 
HI data on inclusive diffractive cross-section, described in section 5.3 of 



5. The HI 2007 Jets parameterization is a Q^-dependent set based on a tune to HI data. 
This fit uses measurements of both the difractive dijet cross-section presented in [27] and the 
inclusive diffractive cross-section presented in 



PDF sets 3,4 and 5 above are next to leading order (NLO) sets, which may make them less 
suited for MC applications. A leading order (LO) gluon might be more stable at small z when 
evolving to lower scales and a LO set will attach better to the LO matrix elements of PYTHIA. 
A LO fit from the HI collaboration [28] is due to be added to the list of PDF sets. 



Parton distributions, by default, are normalized so that they obey the momentum sum rule 



1 

zj{z) dz = 1. 





The motivation for this to hold is described in [2S]. However, since the Pomeron is not a 
physical particle, DPDFs do not implement momentum sum rules. Those from HI add up to a 
momentum sum of roughly 50%. Only the product of the Pomeron flux and the Pomeron PDF 
is meaningful, allowing arbitary separate normalizations of the Pomeron flux and the Pomeron 
PDF. HI choose to normalize their flux so that 

fp/p{xp,t) = 1 when xp = 0.003. 

The standard perturbative multiple interactions framework for pp collisions provides parton- 
parton interaction cross-sections at a fixed CM energy. To turn these cross-sections into prob- 
abilities, one needs an ansatz for the Pomeron-proton total cross section. The single diffractive 
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cross-section is given by equation [131 

oSd = / dxpdt fp/p{xp,t) app(M^ = xps) . gap survival (13) 

(Tpp (effective) 

The normalization of the Pomeron flux {fp/p{xp, t)) is arbitrary and ctsd is parameterized from 
Regge theory. Then crpp(effective) is adjusted accordingly. The value of o"pp often quoted in 
literature is around 2mb |21]. In PYTHIA the default value of cipp (effective) is 10 mb, which 
takes into account screening effects. This value is also used for multiple interactions in diffrac- 
tive systems as described below. It is the main free tunable parameter in high-mass diffraction, 
and along with the choice of Pomeron PDF, can be fitted to represent diffractive event-shape 
data such as average charged multiplicity. The gap survival probability depends on the energy 
of the collision. The higher the energy, the greater the probability of multiple interactions in 
the same event that suppress the rapidity gap. 

To describe the dampening of the perturbative jet cross-section at px -^ by colour screening, 
the actual cross-section is multiplied by a regularization factor Pt/{Pto +Pt)^- Pto is a free, 
tunable parameter of the order 2-4 GeV. The energy dependence of pro is given by 

/ AT \ ecmPow 

/ ecmiNlow \ 



PToiecmm^) = pARei) l-^-^ 



where "ecmNow" is the current energy scale, "ecmRef ' is an arbitrary reference energy at which 
PT(Ref) = Pto (ecmRef) is defined and "ecmPow" is the energy rescaling pace. 

Integrating equation [H] gives the total minijet pair cross-section. The average number of 
jet pairs in an event is given by total mimjet^pair^^s-scction ^ Therefore, increasing the value of 
(Tpp (effective) will reduce the multiple interactions activity per event. This also explains the 
choice of crpp (effective) above. 

Even at a fixed CM energy, the diffractive (high) mass spectrum can lie in the range 10 GeV < 
M < ^/s, with a varying set of parameters (such as the pt cut-off parameter (pro)) along the 
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range. Therefore, multiple interactions are initialised for a few (currently five) different diffrac- 
tive mass values across the range, and all relevant parameters are interpolated between them 
to obtain the behaviour at a specific diffractive mass. Additionally, AB — > XB and AB — t- AX 
are initialized separately. This allows for different beams (or PDFs) on both sides. This also 
facilitates double diffraction. 



2.4 PYTHIA 8.130 vs PYTHIA 6.214 and PHOJET 1.12 

A study comparing the pseudorapidity {r]), transverse momentum (pr) and charged particle 
density {dN^i^/dT]) distributions in PYTHIA 8.130, PYTHIA 6.214 and PHOJET 1.12 at CM 
energy 7 TeV is shown below. Only the SD spectra are compared. 
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Figure 12: ry distribution for SD events at 7 TeV using PYTHIAS. 

It is clearly seen that the addition of hard diffraction to PYTHIA shows an improvement in 
the pt and multiplicity tails, giving a description comparable to PHOJET, which also has hard 
diffractive scattering. 
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Figure 13: pr distribution for SD events at 7TeV using PYTHIAS. 
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Figure 14: multiplicity distribution for SD events at 7TeV using PYTHIAS. 
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3 Conclusions and future outlook 

Diffraction is not well understood, and the method employed in describing diffractive processes 
in PYTHIA is only one among several that have been proposed. This approach to sub-dividing 
the Pomeron-specific parts of the generation into independent sections may not necessarily 
represent the sequence of events in reality. It is important to view the effects as a convolu- 
tion of factors. For example, the total diffractive cross-section is the effect of convoluting the 
Pomeron flux with a Pomeron-proton total cross-section. Neither the Pomeron flux, nor the 
total Pomeron-proton cross-section are known from first principles. This leads to a significant 
uncertainty in the flux factor. 

The value of the assumed Pomeron-proton effective cross-section used in PYTHIA is 10 mb. 
Increasing this value reduces the multiple interaction activity per event but if increased too 
much, pto will be adjusted downwards to ensure that the integrated perturbative cross-section 
stays above the assumed total cross-section. This is the main tuneable parameter in high-mass 
diffraction. 

To further complicate this picture, it is possible that an event that involves a Pomeron-proton 
collision that could have given a diffractive event, in addition, also involves normal multiple 
interactions. This would lead to a topology without rapidity gaps ISU]- Experimentally, such 
events are not triggered as diffractive events. 

A point worth mentioning is that in PYTHIA only the Schuler-Sjostrand description of the 
Pomeron flux includes a separate behaviour of t distribution for double diffraction. Since dou- 
ble diffractive events do not have an outgoing proton, it is experimentally difficult to measure 
the t distribution. 

Central diffraction, although tiny, contributes to the total pp cross- sect ion. Its addition to 
PYTHIA can be foreseen in the not so distant future. 
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With the inclusion of this perturbative description of diffraction in PYTHIA, the kinematic 
predictions of PYTHIA in the diffractive areas is now comparable with PHOJET, which uses 
a different but related physics model called the Dual Parton Model (DPM) [31]. Although the 
diffractive part of PYTHIA has made considerable progress, there are still some issues that need 
addressing. Most importantly, comparisions with data have to be made in order to validate the 
model. 
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